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Abstract 

We present the first evidence from lattice simulations that the magnetic monopoles in 
three dimensional compact quantum electrodynamics (cQEDs) with Nf = 2 and Nf = 4 
four-component fermion flavors are in a plasma phase. The evidence is based mainly 
on the divergence of the monopole susceptibility (polarizability) with the lattice size at 
weak gauge couplings. A weak four-Fermi term added to the cQEDa action enabled 
simulations with massless fermions. The exact chiral symmetry of the interaction terms 
forbids symmetry breaking lattice discretization counter-terms to appear in the theory's 
effective action. It is also shown that the scenario of a monopole plasma does not depend 
on the strength of the four-Fermi coupling. Other observables such as the densities of 
"isolated" dipoles and monopoles and the so-called specific heat show that a crossover 
from a dense monopole plasma to a dilute monopole gas occurs at strong couplings. The 
implications of our results on the stability of U (1) spin liquids in two spatial dimensions 
are also discussed. 
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1 Introduction 



Gauge field theories play an important role in both high energy and condensed mat- 
ter physics. The mechanism of quark confinement in gauge theories with dynamical 
fermions such as QCD remains one of the most elusive subjects in particle physics. As 
a result, model field theories play a significant role in studying this phenomenon. Three 
dimensional parity-invariant compact quantum electrodynamics is such an interesting 
and challenging field theory with rich dynamics that resemble four-dimensional QCD. 
It is an asymptotically free theory, because the gauge coupling has mass dimension 
one and thus provides the theory with a natural scale that plays the role of Aqcd in 
four dimensions. Polyakov in his pioneering work on quenched CQED3 |1] showed an- 
alytically that static electric charges are confined via a linear potential for arbitrarily 
small values of the gauge coupling. More specifically, he showed via duality of electric 
and magnetic monopole-like instanton charges that the model is equivalent to a three- 
dimensional Coulomb monopole gas described by a sine-Gordon effective action; this in 
turn leads to a nonzero photon mass and area law for the Wilson loop [2]. 

The situation is less clear when CQED3 is coupled to Nf massless four- component 
fermionic flavors, because the interaction between monopoles and antimonopoles is 
changed by the vacuum polarization. A simple way of seeing why massless fermions 
might be expected to have a dramatic effect is to observe that as a result of the Dirac 
quantization condition the combination eg {g is the charge of the magnetic monopole) 
is a renormalization group (RG) invariant [3J. Given that the renormalized electron 
charge en < e due to screening by e~e+ virtual pairs then the renormalized monopole 
charge gn > g. Hence, virtual e~e^ pairs antiscreen the monopole-antimonopole {mm) 
interaction. If the monopoles are in plasma phase at least for small Nf values, then 
based on the dual superconductor model [U |5] the electric charges are linearly confined. 
In gauge field theories the particles of the vacuum that are analogous to the electrons 
of the superconductor are the magnetic monopoles. The monopoles set up magnetic 
currents which confine the electric field between the charges into a narrow flux tube, in 
a similar way to the electric currents around magnetic flux tubes in an ordinary super- 
conductor. Since this narrow flux tube has a constant energy/length, it gives rise to a 
linearly confining potential. 

The issue of the (non-)existence of a monopole plasma phase in CQED3 coupled 
to Nf massless fermionic flavors has been addressed analytically by various authors. 
Different approaches often based on perturbative renormalization group (RG) analysis 
of an approximate dual anomalous sine-Gordon (ASG) action led to different results 
depending on the type of approximations in the calculations. Using an electrostatic 
argument and an RG calculation the authors of [6] claimed that the interactions among 
magnetic dipoles screen the logarithmic mm potential for arbitrarily large but finite Nf 
back into the Coulomb form at large distances. This result was confirmed by a self- 
consistent variational analysis of the dual ASG theory [?]. The results of [SI [7] were 
criticised by the authors of [8] who showed in a systematic RG analysis that for large 
Nf the monopole operators are irrelevant in the infrared limit and the physics of the 
system is controlled by a conformally invariant fixed point (in the context of cuprate 
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superconductors discussed later in this section it is known as the algebraic spin liquid). 
Arguments based on analysis of topological symmetries [9] produced results consistent 
with ^ . In addition, the authors of [10] claimed that for Nf > 2 the average size of the 
mm dipoles collapses to zero leading to non-compact QED3, provided the fermions are 
massless. If the fermions have a small mass then the monopoles are in a dipolar phase. 
In a more recent RG calculation the authors of [11] claimed that for Nj > A^^'''* = 36 the 
fermions are deconfined, for 20 < Nf < 36 they can be either confined or deconfined, 
depending on the monopole density and for Nf < 20 the fermions are confined. 

Lattice simulations provide a reliable non-perturbative tool for studying the role of 
magnetic monopoles in CQED3. So far, there have not been any simulations that address 
directly the (non-) existence of a monopole plasma phase in CQED3. The inclusion of 
massless fermions in the compact U{1) gauge action makes simulations difficult due to 
the non-local interactions generated when integrating over fermionic variables. There- 
fore, the authors of [12j addressed the issue of electric charge confinement in CQED3 via 
lattice simulations of an effective U{1) lattice gauge theory with a variety of nonlocal 
interactions in the time-like direction that mimic the effects of gapless/gapful matter 
fields. The main result of [12] is that for certain power-law decaying interactions (mim- 
icking coupling to massless matter fields) a second order phase transition separates a 
confined from a deconfined phase. The existence of a deconfined phase in the effective 
theory indicates that when CQED3 is coupled to a large number of massless matter fields 
the theory may be in the deconfined phase. It has been also shown with Monte Carlo 
simulations [13] that charged particles with ln(r) interactions exhibit a phase transition 
at a critical temperature between a dilute dipole gas and a monopole plasma. This 
result also provides indirect evidence that monopoles in CQED3 may be in a dipolar 
phase above a certain A^^'''*. 

In this paper we present the first attempt to resolve the controversy in the analytical 
literature via lattice simulations of cQED3 with Nf = 2 and Nf = 4. Massless fermion 
simulations were enabled with the inclusion of a weak (unable to break chiral symmetry 
on its own) four- Fermi term in the theory's action. The results that are largely based on 
the diverging monopole susceptibility (polarizability) with the lattice extent L at weak 
couplings imply that the monopoles are in a plasma phase. The details of the lattice 
model including the role of the four- Fermi interaction are discussed in Sec. 2. Recent 
simulations of non-compact QED3 (ncQED3) with an extra weak four-fermi term [H] 
showed that the magnetic charges (which unlike in CQED3 they are not classical solutions 
of the theory) form tightly bound mm dipoles, because in this case the Dirac strings 
carry a non-vanishing contribution to the pure gauge (non-compact) part of the action 
|15j . The different monopole dynamics in lattice CQED3 and ncQED3 for A^^ < 4 at 
weak gauge couplings imply that the two models have different continuum limits. The 
authors of [161 CI] performed simulations of both CQED3 and ncQED3 with Nf = 2 at 
strong gauge couphngs and concluded that the two formulations may be equivalent. This 
suggestion was largely based on comparisons of chirally extrapolated data for the chiral 
condensate versus monopole density which appeared to collapse on the same curve for 
the two QED3 formulations. A similar claim was presented for weak couplings [T7] based 
on simulations with a lattice size L = 32. These results, however, are questionable given 
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that Nf = 2 ncQED3 simulations with L = 50 [H] and later on with L = 80 |T9] did 
not provide any evidence for the existence of a nonzero chiral condensate. The principal 
obstruction to a definite answer in ncQEDs is the large separation of scales in the theory, 
i.e. the fermion dynamical mass is at least an order of magnitude smaller then the natural 
cutoff scale e^. In addition, large finite volume effects resulting from the presence of a 
massless photon in the spectrum prevent a reliable extrapolation to the thermodynamic 
limit. So far, the evidence from lattice simulations of ncQEDa is that A^™* < 1.5 
fl9\ |20] . It should be noted though, that recent lattice simulations of ncQEDa with an 
additional weak four-Fermi term [H] hinted at evidence that chiral symmetry may be 
broken up to Nf = 4. Analytical results based on self-consistent solutions of Schwinger- 
Dyson equations (SDE) [211 |22] claimed that to detect chiral symmetry breaking for 
Nf > 1.5 lattice volumes much bigger than the ones currently used in simulations are 
required. The most recent SDE-based analytical calculations showed that Nf^^'^ ^ 4 
f2T\ 123] • In addition, a gauge invariant calculation based on the divergence of the chiral 
susceptibility resulted in A^^"* = 2.2 ^3] and an RG approach on QED3 with extra 
irrelevant four- Fermi interactions resulted in N™^ = 6 [25j. An alternative approach 
based on simulations of the (2 + l)d Thirring model at infinite coupling which may 
belong to the same universality class as ncQEDs resulted in Nj^^^ = 6.6(1) |26] . 

Although cQEDa is a model field theory used for studying elementary particle physics 
phenomena, and is also an interesing basic field theory on it own right, it acquires more 
concrete phenomenological significance in condensed matter physics, because (2 + 1)0? 
field theories with a compact U{1) gauge field coupled to gapless relativistic fermions 
arise as low energy effective field theories in two dimensional strongly correlated electron 
systems, such as the cuprate superconductors [271 1211 |29l [30]. Strong interactions lead 
to correlated electron motion resulting in unconventional states of matter with "frac- 
tionalized" quantum numbers where the quasiparticle approach of Landau's Fermi liquid 
theory is not valid. It is well-known from experiments that, for cuprates, the Mott in- 
sulating state at zero doping is the Neel antiferromagnetic state but the nature of the 
connection between the undoped Mott state and the doped rf-wave superconductor is still 
under theoretical debate. In a specific incarnation of Anderson's [31] resonating valence 
bond idea, it was proposed [2S1 EH] that the so-called U{1) spin liquids are the phases 
of matter that play a significant role in understanding underdoped cuprates. These fea- 
tureless quantum paramagnetic states with no broken symmetries or long-range order, 
also known as critical or algebraic spin liquids (ASL), behave as if the system is at a 
critical point without the fine-tuning of any parameter. The physical picture can be 
visualized in terms of valence bonds between pairs of spin singlets separated by arbi- 
trarily large distances and the unpaired charge- neutral gapless spin 1/2 spinous interact 
strongly with the valence bond background. The flux of the emergent U{1) gauge field 
arises from extra topological conservation laws not present in the microscopic theory 
[32] . In a low energy description of U{1) spin liquids, the spinous with a linear dis- 
persion are coupled minimally to an emergent compact U{1) gauge field, resulting in 
cQED3. Of particular importance to the physics of cuprates is the so-called staggered 
flux (or (i-wave resonating valence bond) phase which is formally described by CQED3 
(the anisotropy of the interactions is neglected) with Nf = 2 four-component fermions. 
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Another ASL with an important role in strongly correlated electron systems is the so- 
called TT-flux state, which is described by Nf = A CQED3. Ghaemi and Senthil |33] showed 
that the staggered flux spin liquid state may be connected to the Neel antiferromagnetic 
state via a second order quantum phase transition. For extensive reviews on high Tc 
superconductivity resulting from doping a Mott insulator we refer the reader to IlM] and 
references therein. These spin liquid phases are stable provided the magnetic monopoles 
are not in a plasma phase and hence are unable to induce linear spinon confinement. In 
addition to the works mentioned in a previous paragraph regarding the role of magnetic 
monopoles in cQEDs [HI [71 El El [TOl [H] , Alicea |35j found a monopole operator which 
represents a symmetry allowed perturbation, and speculated that this may destabihze 
the staggered flux phase leading to charge confinement. The numerical evidence pre- 
sented in this paper in favor of the existence of a monopole plasma phase for Nj < 4 
implies that both the staggered-fiux and the vr-fiux spin liquids are unstable to spinon 
linear confinement. This in turn leads to Neel antiferromagnetic order where the chiral 
condensate (i'lp) corresponds to the staggered magnetization. It should also be noted 
that in a more phenomenological approach an anisotropic version of ncQEDa has been 
proposed for the phase fluctuations in the pseudogap phase (also known as algebraic 
Fermi liquid phase) of cuprate superconductors [36] . 

This paper is organized as follows: In Sec. 2 we present the lattice model, the various 
monopole observables, and the simulations parameters. In Sec. 3 we present and discuss 
the simulation results and in Sec. 4 we present our conclusions and also point to possible 
future expansions of this project. 



2 The Model 

In this first numerical exploration of the role of magnetic monopoles in CQED3 with 
an additional four-Fermi term we have chosen the simplest Z2 chirally symmetric four- 
Fermi interaction which for practical purposes is preferable over terms with a continuous 
chiral symmetry, because the latter are not as efficiently simulated due to the presence of 
massless modes in the strongly cut-off theory. For computational purposes it is useful to 
introduce the auxiliary field a = glxiXi (summation over the index i is implied), where 
Xi, Xi are ^-component staggered fermion fields and gl is the four- Fermi interaction 
coupling. The lattice action in terms of real-valued link potentials 6^1, is given by the 
following equations: 

1 = 1 ^ x,x' ' X X,fJ,<U 

where 

e^,, = e^,{x) + 0,{x + fi)-e^{x + u)-0,{x), (2) 

Q{x,x') ^ (3) 

fJ. (x,x) 
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The indices x, x' consist of three integers (xi, X2, X3) labelhng the lattice sites, where the 
third direction is considered timehke. The symbol {x, x) denotes the set of the eight dual 
lattice sites x surrounding the direct lattice site x, where the a field lives [37]. The r]^{x) 
are the Kawamoto-Smit staggered fermion phases (—1)^1"' '"^a'-i^ designed to ensure 
relativistic covariance of the Dirac equation in the continuum limit. The inverse gauge 
and four- Fermi couplings are given by /3 = ^ and f3s = respectively, and a is the 
physical lattice spacing. The boundary conditions for the fermion fields are antiperiodic 
in the timelike direction and periodic in the spatial directions. In the weakly coupled 
(/3 — )■ 00) long- wavelength limit eq. ([T]) describes Nf four- component Dirac fermions |38j . 

Performing simulations with massless fermions even with the reduced Z2 chiral sym- 
metry has substantial advantages, both theoretical and practical. The theory has the 
exact symmetry of the interaction terms, which forbid chiral symmetry breaking coun- 
terterms from appearing in its effective action. In addition, because of the large nonzero 
vacuum expectation value of the a field at strong gauge coupling^ or its fluctuations 
at weak couplings, the Dirac operator is nonsingular even with m = and its inversion 
is very fast. Another advantage of simulations with m = 0, is that we do not have 
to rely on often uncontrolled m — )■ extrapolations to calculate various observables in 
chiral limit. For these reasons both the non-compact and compact lattice versions of 
the theory have been successfully simulated in (3 -|- l)d [391 HO] and showed that QED4 
is a logarithmically trivial theory and the systematics of the logarithms follow those of 
the Nambu— Jona-Lasinio model rather than those of the scalar as often assumed. 
Unlike (3 + l)d where the four- Fermi term is a marginally irrelevant operator, in {2 + l)d 
it is a relevant operator. It is well-known that the (2 -|- l)d Gross-Neveu model (GNM3) 
although non-renormalizable in the weak coupling perturbation theory, is renormalizable 
in the 1/Nf expansion ^T]. At sufficiently strong couplings chiral symmetry is sponta- 
neously broken, leading to a dynamically generated fermion mass S ^ (a) » m. The 
interacting continuum limit of the theory may be taken at the critical coupling g"^^ (at 
which the gap S/A[/y — )■ 0), which defines an ultraviolet-stable renormalization group 
fixed point. 

In QED3, as the gauge coupling is varied (with the four- Fermi coupling kept fixed 
at some weak value (7^ < g^J, depending on the value of Nf the model is expected 
to undergo either a chiral phase transition or a sharp crossover from a strong coupling 
phase (where (xx) 7^ 0) to a weak coupling phase where (xx) is either zero or very small 
and possibly undetectable in lattice simulations. Hereafter, we will use the term "chiral 
transition" to denote either a chiral phase transition or a sharp crossover from strong to 
weak gauge couplings. Near the transition the weak four-Fermi term is expected to play 
a dominant role as compared to the ultraviolet-finite gauge interaction. Simulations 
of ncQEDa [14J showed that the order parameter scales with critical exponents close 
to those of GNM3 and the scaling region is suppressed by a factor ~ gg. The GNM3 
scaling is expected to be valid for both compact and non-compact lattice formulations. 
It should also be noted that in the large- A^/ and /3 — 00 limits the four-Fermi term is 

^At strong couplings, pure QED3 simulations are dramatically slowed down by the strong gauge field 
fluctuations. 
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an irrelevant operator in the RG sense [25]. 

The simulations were performed with the standard Hybrid Molecular Dynamics 
(HMD) R algorithm. We used conservatively small values for the HMD trajectory time- 
step dt and ensured that any 0{dt^) systematic errors are smaller than the statistical 
errors for different observables. For lattice sizes smaller than or equal to 24^ we used 
dt = 0.005 and an HMD trajectory length r = 1 and for 32'^ we used dt = 0.0025 with 
a trajectory length r = 2. 

The magnetic monopoles in the lattice model are identified following the standard 
DeGrand and Toussaint approach [12]. The plaquette angles 6^,^ are written as 

e^,^ = ef,u + 2'Ks^^{x), (4) 

where Q^i, lie in the range (— vr, tt] and s^^{x) is an integer that determines the flux due 
to a Dirac string passing through a plaquette. The gauge invariant integer number of 
monopole charges on the dual lattice sites x are then given by 

M(£) = e^,x^+s,x{x), (5) 

where is the forward lattice derivative and M G {0, ±1, ±2}. Since on a three-torus 
the total magnetic charge M[x) = we define the density of monopole charges as 

PM = ^X^|M(x)|. (6) 

The mere counting of monopoles does not provide any useful information on whether 
their presence has any impact on the model's confining properties. As already discussed 
in Sec. 1 a monopole plasma is required for linear confinement of electric charges. 

The observable that provides information on whether the monopoles are in a plasma 
or a dipolar phase is the monopole susceptibility Xm [43] : 

X™ = -$^(r'M(0)M(r)). (7) 

r 

The susceptibility is the polarizability of the monopole configurations; this can be readily 
seen by adding a uniform magnetic field term —B rM{r) to the dual monopole action 
[l3] and evaluating Xm = hi Z / B\b=oi where Z is the monopole partition function. 
If the magnetic charges are in a plasma phase, then Xm diverges with the lattice size L, 
implying that external magnetic fields are shielded. A finite Xm means that monopoles 
and antimonopoles form a polarized gas of mffi dipoles, which is what was observed 
in ncQEDa simulations [T?]. The situation may be very different in CQED3, where the 
monopoles are classical solutions of the theory and they may exist in a plasma phase at 
least for small Nf values. Results from numerical simulations of CQED3 with Nf = 2 
and Nf = 4 presented in Sec. 3 favor the existence of a monopole plasma phase. Also, 
as shown in |11] , in the infinite volume limit further manipulations lead to a form of Xm 
expressed as a Fourier transform of a two-point correlation function at zero wavevector: 

Xm oc ^(0^^(x)e^^(O)) oc ^(s^^(x)s^^(O)). (8) 
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The observable Xm has been rarely measured in simulations with dynamical fermions, 
because it is very noisy due to near cancellations of monopole-monopole and monopole- 
antimonopole contributions. With the inclusion of the four-fermi term in the QED3 
action the algorithm became very efficient and Xm has been measured with an acceptable 
signal-to-noise ratio at weak gauge couplings. We generated ~ 10^— 2 x 10^ configurations 
for the largest L = 32 lattice and ^ 3 x 10^ — 7 x 10^ configurations for the smaller 
lattices (L = 8, 24). 

We also measured Xi given by 

Xi = -(M(0)M(1)), (9) 

which includes the contributions in Xm from adjacent lattice cubes only. In addition, we 
measured X2 given by 

X2 = - $^(r2M(0)M(r)), (10) 

which includes the terms of Xm where two neighboring magnetic charges share either a 
cube face, an edge or a corner. Xi ~ X2 ~ Xm indicates that the main contribution to 
Xm comes from tightly bound mfh pairs. Therefore, a comparison of Xm with xi and X2 
provides information on whether the monopole configurations are dominated by tightly 
bound dipoles or not. This will become clearer in Sec. 3 where we compare the behavior 
of these observables (as a function of L) for both cQEDs and ncQEDa. 

3 Results 

In this section we present results from simulations of the lattice model in eq. ([1]) with 
Nf = 2 and Nf = 4 fermion fiavors. Before presenting data for the monopole observables 
we present results for the chiral condensate (xx) versus /3 near the Nf = 4 strong 
coupling chiral transition. In the infinite gauge coupling limit (/3 — )■ 0), it is known 
rigorously that chiral symmetry is broken [ISj for values of Nf below a certain critical 
value. Simulations of QED3 with staggered fermions and /3 = showed that the theory 
undergoes a second-order phase transition at Nf ^ 8 [?6]. With the extra weak four- 
Fermi term the infinite gauge coupling transition is shifted towards larger Nf values, 
depending on the value of Ps- Therefore, in cQEDs as P increases and Nf is larger 
than a putative A^j"* (for Nf > Nj"^ the monopoles are in the dipolar phase) there 
must exist a chiral phase transition at some critical value of the gauge coupling Pc- For 
Nf < A^™'*, since the chiral order parameter is small at weak gauge couplings, the relic 
of the transition may persist as a sharp crossover between weak and strong couplings 
with a tail of exponentially suppressed xX extending to weak couplings. The bulk of 
the simulations presented in this paper were performed with a fixed Ps = 2.0, which 
is larger than the critical coupling = 0.835(1) for the three-dimensional Nf = 4 
GNM3 (and hence larger than P"^^ for Nf = 2 GNM3). Therefore, the four- Fermi term 
with Ps = 2.0 cannot break chiral symmetry on its own. However, as already mentioned 
in Sec. 2, in three dimensions the weak four- Fermi term is expected to play a dominant 
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Figure 1: (Color online) (|xxl) vs /3 extracted from simulations with Nj = 4 on 24^ and 32'^ lattices. 
The sohd curve represents the fitting function (eq. (HH)) for f3 E [1.025, 1.200] . 

Table 1: Values of (|xx|) for Nf = 4 from simulations on 24"^ and 32^ lattices near the crossover. 





(|XX|)24 


{\XX\)32 


1.30 


0.0612(9) 


0.0658(17) 


1.35 


0.0372(11) 


0.0383(15) 


1.40 


0.0236(3) 


0.0238(10) 


1.45 


0.0178(5) 


0.0154(9) 


1.50 


0.0134(5) 


0.0125(9) 


1.60 


0.0112(2) 


0.0080(2) 



role near the chiral transition as compared to the uhraviolet-finite gauge interaction. It 
was shown in simulations of ncQEDa |T4] that the critical exponents extracted for the 
Nf = 4 transition are close to the GNM3 ones, and small deviations hinted at evidence 
for nonzero fermion mass generated by the gauge field dynamics. On a finite volume 
lattice near the transition the values of XX may change sign due to tunnelling events 
between the Z2 vacua resulting (xx) = 0. In order to take into account these tunnelling 
events and following similar analyses of the Ising model we measured the effective order 
parameter instead of (xx)- We fitted the data extracted from simulations with 

L = 24 for (I XX I) versus /3 to the standard scaling relation for a second-order phase 
transition order parameter 

{\XX\) = A{P,-Pf-. (11) 

For the fitting range /3 G [1.025, 1.200] we obtained Pc = 1.35(1) and Pm = 0.94(3) 
with an acceptable fit quality given by x^/DOF = 1.8. The extracted value of (3m is in 
very good agreement with the (3^ = 0.93(3) of iV/ = 4 GNM3 |17]. The GNM3 scaling 
confirms our earlier assertion that the exact chiral symmetry of the lattice action forbids 
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Figure 2: (Color online) Binder's cumulant vs (3 ior L — 16, 24, 32 near the Nf = 4 chiral crossover. 



symmetry breaking lattice discretization counterterms from appearing in the model's 
free energy. The 24^ data (together with some 32^ data) and the fitted curve are shown 
in Fig. [H By expanding the fitting window towards larger values of /3, the value of Pm 
increased and the fit quality deteriorated. The fit quality deteriorated dramatically when 
data points above (3 = 1.35 were included. If, however, (3 = 1.35 were a critical coupling 
then one would expect the effective order parameter to obey the finite size scaling relation 
(Ixxl) ~ L'^"'^" - L-°-^^ {/3m/iy = 0.927(15) in Nj = 4 GNMg gTj). However, as shown 
in Table 1, instead of observing a decrease of with L at the putative critical 

coupling 13 = 1.35 we observe that for f3 = 1.35, 1.40 the values of the effective order 
parameter on 24^ and 32^ are equal within statistical errors. The values of for 
P = 1.45, 1.50 from the two lattices also agree within 1-2 standard deviations. Finally 
at /3 = 1.60 which corresponds to the weakest gauge coupling in Table 1 large finite 
volume effects (the physical volume shrinks with (3) result in (|xxl)32 < (|xxl)24- These 
observations together with the failure of eq. f lTT]) to provide acceptable fits when weak 
couplings were included in the fitting window constitute serious evidence that instead of 
a chiral phase transition we have a crossover from strong to weak couplings. It should 
be noted that in ncQEDg although we observed small deviations of the exponents from 
the GNM3 ones, above the transition there was a decrease of with L, possibly 

because any tiny non-zero chiral condensate is "swallowed" by finite volume effects. This 
difference between the two lattice QED3 formulations could be attributed to the different 
magnetic monopole dynamics, i.e. in CQED3 the monopoles may be in a plasma phase, 
which in turn leads to an enhanced chiral condensate. Before concentrating on the role 
of monopoles in cQEDs, we will study the behavior of the Binder cumulant Ub{f3, L) [18] 
defined by 

, iixxr) 



(12) 
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Figure 3: (Color online) Xm vs. L for (a) Nf — 2, and (b) Nf = 4. The fitting functions (eq. (|13l) ) are 

represented by solid lines. 



and measured on different lattice sizes near the chiral crossover. Near a second order 
phase transition and with sufficiently large lattices (where subleading corrections from 
finite L are negligible) f/f, = fb{{Pc — /3)L^^'^)- Therefore, at a critical coupling the lines 
connecting data of the same L are expected to cross at a universal value Ub = U^. In 
a symmetric phase as L ^ oo t/5 ^ 0. For Nj = A GNM3 = 0.232(8) In 
Fig. |2]we present Ubi/S, L) data for Nf = A cQEDg with L = 16, 24, 32. Although Ub is a 
lot noisier than the effective order parameter it is clear that the lines joining data with 
the same L do not cross for /3 < 1.45. This observation provides additional evidence in 
favor of the crossover scenario instead of a phase transition for /3 < 1.5. The crossings 
of the constant L Ub lines occur at /3 ~ 1.6 where the values of Ub are less than 0.02. As 
stated earlier even if for /3 > 1.6 chiral symmetry is broken, the physical lattice volume 
is smaller and therefore it is plausible that finite size effects make the phase look as if it 
is symmetric. 

Next, we turn our attention to the role of magnetic monopoles. As discussed in Sec. 2, 
the most relevant observable for deciding whether the monopoles are in a plasma or a 
dipolar phase is the monopole susceptibility Xm- Its short distance contributions xi and 
X2 elucidate further the situation. It is instructive to compare data for these observables 
from both the compact and non-compact lattice formulations of QED3. Recent ncQEDa 
simulations [H] showed that the magnetic charges form tightly bound mm dipoles. In 
Figures [3] (a) and (b) we present the Xm results at weak gauge couplings for Nf = 2 and 
Nf = 4, respectively. The data are fitted to a power-law relation 

Xm = cL- (13) 
and the extracted values of the exponent a are shown in Fig. |H The clear increase of 
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Figure 4: (Color online) Exponent a vs. /3 extracted from fits with eq. (fT3|) for Nf — 2 and Nf — 4. 



Xm with L in Fig. [3] and the non-zero values of the exponent a in Fig. H] imply that the 
magnetic charges are in a plasma phase. As expected, the values of a for Nf = 2 are 
larger than those for Nf = 4. This can be understood in terms of the renormalization 
group invariant Dirac quantization condition eg = enQR = n/2 (n is an integer and 
the subscript R denotes renormalized charges): As Nf increases the e^e~ interaction 
decreases due to enhanced screening from virtual fermion-antifermion pairs, which in 
turn implies that the mfh interaction is antiscreened. In addition the Nf = 4 values 
of a appear to decrease with (3. It is possible that for large Nf values, the extracted 
values of a may be affected by finite size effects. This can be understood as follows: The 
interaction among magnetic dipoles leads to a screening of the mm interaction |6l [7], 
resulting in a small density of unbound magnetic charges. This mechanism implies that 
the existence of a monopole plasma may be a very long distance effect. Therefore, 
simulations on larger lattices may be required in order to extract more accurate values 
of a. The existence of the monopole plasma depends solely on the gauge field dynamics 
and it is not expected to depend on the four-Fermi coupling provided the latter is 
weak enough. This is supported by the f3s > 2.0 data in Fig. [3] (b). The data for 
Ps = 8-0, 12.0, 16.0 from simulations with L = 8, 24 almost collapse on a single curve 
and the values of a for the three different (3s are consistent with a = 0.61(7) which is 
close to the (3s = 2.0 value a = 0.73(3). The slightly smaller value of a in the (3s ^ oo 
limit could be attributed to the fact that the additional four-Fermi term enhances the 
e~e~^ interaction, implying an enhanced antiscreening effect on the mm interaction. 

At this point it worth comparing data for Xm from both CQED3 and ncQEDs. In 
Fig. 13 we present Xm versus L for ncQEDa above the chiral transition/crossover for 
Nf = 2,4, (3s = 2.0 and (3 = 0.4,0.5. The horizontal lines give excellent fit qualities, 
implying that the values of Xm do not depend on L, because the magnetic charges are 
in the dipolar phase. 

In Figures [6] (a) and (b) we plot Xm together with its short distance contributions Xi 
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Figure 5: (Color online) Xm vs. L for ncQEDa with Nf = 2,4 and /3 = 0.4, 0.5. The horizontal lines 

represent fits to the data. 




Figure 6: (Color online) Xm, Xi, and X2 vs. L for (a) cQEDa with Nf = 2,(3 ■ 

and (b) ncQEDg with Nf = 2,4 and (3 = 0.4. 
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Figure 7: (Color online) /?^(xx) /3m for A^y = 2, iV/ = 4 at different values of the coupling /?. The 
physical volume {L/ (3)^ is constant for each Nf value. 



and X2 versus L for CQED3 and ncQEDs, respectively. For CQED3, the Xi and X2 data fall 
on different horizontal lines below Xm-, because the divergence of Xm which is a signature 
for a monopole plasma comes from long distance contributions. In contrast to this, the 
ncQEDs data for xii X2 and Xm coincide within statistical errors. This confirms that 
in ncQEDa the contribution to the polarizability comes from tightly bound mm dipoles 
which in the continuum (/3 — 00) limit may disappear by collapsing into zero size. 

Next, we check whether the gauge couplings /3 = 2.0 and /3 = 1.8 are in the asymp- 
totic scaling regimes for Nf = 2 and Nf = 4, respectively. Since for an asymptotically- 
free field theory the ultraviolet behavior is governed by the gaussian fixed point at the 
origin, then the continuum limit of the model lies in the limit /3 — )■ 00, and all physical 
quantities should be expressible in terms of the scale set by the dimensionful coupling e^. 
To compare simulation results taken at different couplings (lattice spacings), therefore, 
it is natural to work in terms of dimensionless variables such as (3m, L/P, and /3'^{xx)- 
As the continuum limit is approached, data taken at different /3 should collapse onto 
a single curve when plotted in dimensionless units. To check whether lattice data are 
characteristic of the continuum limit we plot the dimensionless chiral condensate (xx) 
versus the dimensionless fermion bare mass Pm for Nf = 2 (with /3 = 2.0, L = 24 and 
P = 2.4, L = 30) and Nf = A (with /3 = 1.8 , L = 24 and /3 = 2.4, L = 32 ) in Fig. [3 
For both Nf = 2 and Nf = 4, the values of P'^{xx) fo^" the two values of P agree within 
1-2 standard deviations, implying that lattice discretization effects are small. 

In order to complete the picture presented in the previous paragraphs we also mea- 
sured various densities: (i) the density pd of "isolated" tightly bound dipoles, which are 
mm pairs on adjacent cubes, and each charge in the dipole does not share any other 
cube face with with a second opposite charge; (ii) the density pi of positive magnetic 
charges that do not share a cube face with an antimonopole; and (iii) the density p2 of 
"isolated" positive magnetic charges, that not share a cube face, an edge or a corner 
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Figure 8: (Color online) p^, pi, and p2 vs. /3 for (a) Nf — 2, and (b) Nf — 4. 



with an opposite charge. We fitted the data for p^, Pi and p2 versus j3 to an exponential 
function = oi exp(— 02/?). The data and the fitted curves for Nf = 2 and Nf = A 
are shown in Figures [8] (a) and(b), respectively. The results are from simulations on 32^ 
lattices and a comparison with data on 24^ lattices showed that finite volume effects are 
negligible. It is clear from the peak values of p2 that a crossover from a dense monopole 
plasma phase at strong couplings to a dilute monopole gas at weak couplings occurs. 
The values of the crossover couplings are /3 ~ 1.25 and /3 ~ 1.10 for Nf = 2 and Nf = 4, 
respectively. The values of 02 are presented in Table [2l The fact that the exponential 
function accurately fits the data implies that there is no phase transition that would 
result in an abrupt decrease of the densities. As expected, the values of the parameter 
02 for the "isolated" monopoles are a bit larger than the respective one in quenched 
cQEDs where 02 ~ 5 [12] , whereas for the "isolated" dipoles the values of 02 are smaller 
than the quenched value a2 = 6.6 |12]. The crossover from a dense monopole plasma to 
a dilute monopole gas is also supported by the behavior of the so-called specific heat Cy 
defined in a way analogous to the specific heat in spin models (with the temperature T 

Table 2: Values of parameter a2 extracted from fits of /(/?) = ai exp(— a2/3) on pd, pi and p2 vs. (3 

data. 
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Figure 9: (Color online) Specific heat C„ vs. /3 for TV/ = 4 and L = 16, 24, 32. 



interchanged with the gauge couphng g"^) as follows: 



,d{P) 



dp 



0' 



\nZ 

a/32 



(14) 



where Z is the lattice partition function and P = ■^'Y^x fj.<vi^ ~ co?,Q^i,x) is the pure 
gauge part of the action per unit volume. In Fig. |9]we plot versus (3 for Nf = 4 and 
L = 16,24,32. It is clear that develops a lattice size-independent peak at /3 = 1.1 
The L-independent peak of implies that a smooth crossover takes place in the gauge 
field dynamics at /3 = 1.1, which coincides with the crossover from a dense monopole 
plasma to a dilute monopole gas. 



4 Conclusions and Outlook 

cQEDs is an interesting field theory due to its similarities with QCD and its close relation 
with QCD-like theories [15]. In this paper we presented the first analysis of monopole 
dynamics in cQEDs with Nf < 4 based on results from lattice simulations. Fast sim- 
ulations with massless fermions were enabled by adding a weak four-Fermi interaction 
to the cQEDs action, because the vacuum expectation value of the a meson field which 
appears explicitly in the semi-bosonized action acts like a fermion mass and makes the 
Dirac matrix inversion fast. In addition, in the presence of the four- Fermi term the action 
has an exact Z2 chiral symmetry, which forbids symmetry breaking lattice discretization 
counterterms in the free energy. 

The monopole susceptibility (polarizability) diverges with the lattice extent, implying 
that the monopoles are in a plasma phase, which in turn leads to linear confinement of 
electric charges. Simulations at a single value of /3 for Nf = 4 showed that the monopole 
plasma scenario does not depend on the four-Fermi coupling Ps when (3^ is sufficiently 
large. The cQEDs results for the monopole susceptibility were contrasted with ncQEDa 
data where Xm is independent of L and its major contribution comes from charges on 
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adjacent lattice cubes, implying that the single lattice spacing size dipoles at finite 
P may collapse to zero size in the continuum limit. In addition, the behavior of the 
effective chiral order parameter for Nf = 4 implies that a crossover instead of a 

transition takes place at strong couplings, which could be an outcome of a l^(r) ~ r 
electric potential. Also, the behavior of the density of "isolated" monopoles favors a 
scenario of a crossover from a dense plasma of monopoles to a dilute monopole gas at 
weak couplings. This scenario is supported by the L-independent peak of the so-called 
specific heat. Our results imply that for Nf < 4 the continuum limits of CQED3 and 
ncQEDs are different, with linear charge confinement for the former and logarithmic 
confinement for the later provided that Nf = 4 is below the ncQEDs N™^. 

Our conclusions are supported by the results of ^0], where it was shown that an 
isolated magnetic charge has an infinite free energy, both for the dynamical and the 
quenched system. Following similar lines, we plan to test the response of the system to 
the insertion of a static dipole. A vanishing free energy gap would confirm our results 
that monopoles are in a neutral plasma phase. 

The existence of a monopole plasma phase in CQED3 has implications in strongly 
correlated electron systems. More specifically U{1) spin liquids (staggered-fiux spin liq- 
uid corresponds to Nf = 2 and vr-fiux spin liquid to Nf = 4) in two spatial dimensions 
which are believed to describe the underdoped Mott insulator regime in cuprate super- 
conductors may be unstable to spinon confinement. It should be noted, however, that 
the anisotropic interactions of these condensed matter systems have been neglected in 
our model. Recent analytical [51] and numerical [52] results of ncQEDs with Fermi and 
gap anisotropics showed that the velocity anisotropy is relevant in the RG sense and its 
increase leads to a decrease of N™^. 

We are currently expanding the cQED3 simulations to larger Nf values. The plan is 
to search for a putative conformally invariant fixed point at A^^'''* where a phase transition 
from a linearly confining phase to conformal deconfined phase may take place. In the 
condensed matter language this critical point would correspond to a deconfined quantum 
critical point [53j, where a phase transition is expected to occur between a phase with 
Neel aniferromagnetic order (at small Nf) and a paramagnetic critical spin liquid phase 
(at large Nf). In addition, as emphasized in [M] (2 + l)d two-color QCD may provide a 
more appropriate description of algebraic spin liquids than CQED3 which only includes 
Gaussian fiuctuations about the mean field and suffers from various limitations in the 
underdoped regime. The study of the phase diagram of (2 + l)d two-color QCD is 
another non-perturbative problem that requires lattice simulations for reaching definitive 
answers. 
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